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Abstract 

We consider the problem of describing the possible spectra of an acous- 
tic operator with a periodic finite-gap density. We construct flows on 
the moduli space of algebraic Riemann surfaces that preserve the pe- 
riods of the corresponding operator. By a suitable extension of the 
phase space, these equations can be written with quadratic irrational- 
ities. 



1 Introduction 

In the study of periodic potentials of the one- dimensional Schrodinger oper- 
ator, we encounter the following problem. A given periodic potential u(x) of 
the Schrodinger operator 

>}i +u{y) (1.1) 



dy 

determines a Riemann surface, called the spectral curve, which is algebraic 
if the potential has a finite number of energy gaps. Conversely, given a hy- 
perelliptic Riemann surface of genus g, there exists a potential u(x) with the 
corresponding spectrum, which can be extended to a solution u(x, t) of the 
KdV hierarchy. However, this potential is not in general periodic. We can- 
not extend the direct spectral transform to all quasi-periodic potentials. For 



"Columbia University, New York, USA, zakharov@math.columbia.edu 



1 



instance, the spectrum of a quasi-periodic potential may have the structure 
of a Cantor set (see [Tj), making it impossible to associate a Riemann surface 
to it. The problem of determining which Riemann surfaces correspond to 
periodic finite-gap potentials did not at first receive much attention. 

The periods of a potential of the Schrodinger operator are expressed as 
certain Abelian integrals on the corresponding spectral curve. The poten- 
tial has a given period if and only if these integrals are integer multiples of 
the period. Therefore, the subset of curves corresponding potentials with 
a given period forms a transcendental submanifold in the moduli space of 
hyperelliptic Riemann surfaces. 

A description of the possible spectra of a periodic Scrodinger operator, 
not necessarily finite-gap, was given by Marchenko and Ostrovsky in [5] in 
terms of the properties of a conformal map of a certain type. However, their 
approach was difficult from a computational viewpoint. An effective solution 
was given by Grinevich and Schmidt in [6], called the method of isoperiodic 
deformations. Using an idea developed in [5], [Sj, the solution consists of 
finding a set of differential equations on the spectral data that do not change 
the periods of the solution. The submanifold corresponding to solutions of a 
given period is then preserved by these flows. By a suitable extension of the 
phase space, the equations of isoperiodic deformation can be written with a 
rational right-hand side, making them convenient for numerical simulation. 

This paper extends this approach to the study of the periodic finite-gap 
densities acoustic operator 



Finite-gap densities of the acoustic operator and the corresponding periodic 
solutions of the Harry Dym equation 



were first constructed by Dmitrieva in [2], [3], pE], using a Hopf transformation 
to the Schrodinger operator introduced in [12]. Periodic solutions of the 
Harry Dym equation were recently shown to be relevant to the Saffman- 
Taylor problem (see [13]), while finite-gap periodic densities of the acoustic 
operator are related to geodesies on the ellipsoid (see [10], [H], [T5]). 

In sections 2 and 3, we recall the spectral theory of the periodic Schrodinger 
operator and the method of isoperiodic deformations. In section 4, we derive 
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a spectral theory for the acoustic operator by relating it to the Schrodinger 
operator, and in section 5 we construct the equations of isoperiodic deforma- 
tion for the acoustic operator. By extending the phase space, we write these 
equations with quadratic irrationalities. The principal result of this paper is 
Theorem 2, which gives the explicit form of the deformation equations. 

2 Periodic finite-gap potentials of the Schrodinger 
operator 

We first recall the spectral theory of the one-dimensional periodic finite-gap 
Schrodinger operator 

L s =-^ + u(y), (2.1) 

where u(y) is a smooth periodic real-valued potential with period II. For 
references, see [TT] . 

We consider two spectral problems for Ls'- 

1. The standard problem in L 2 (IR): 

L s <f = \ip, \ip(y)\ < oo as y -> ±oo, 

2. The Dirichlet problem on a period: 

L s p = Ay?, (p(y ) = (f{y + II) = 0. 

The spectrum of the first problem is continuous and consists of an infi- 
nite number of segments [A , A]], [A 2 , A 3 ], . . ., with A 2 j < A 2 j+i < A 2j+2 - The 
spaces between these segments (A 2j+1 , A 2j+2 ), which may be of zero length, 
are called the energy gaps. The second problem has a purely discrete spec- 
trum dj(yo), with exactly one eigenvalue inside or on the boundary of each 
of the energy gaps, dj(y ) G [A 2 j_i,A 2 j], including the degenerate gaps. 

The principal case of interest is when the potential u(y) has only a finite 
number of energy gaps of non-zero length, such a potential is called finite- 
gap. Let (— oo, Ao), (Ai, A 2 ), . . . , (A 29 _i, A 2ff ) denote the non-trivial energy 
gaps, and let 7^(2/0) denote the eigenvalue of the Dirichlet problem lying in 
the j-th non-trivial gap. 

The direct spectral transform assigns to the finite-gap potential u(y) the 
following data 
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1. A hyperelliptic Riemann surface T of genus g together with a two- 
sheeted covering A : T — > CP 1 ramified at A , . . . , X 2g and oo 

2. A meromorphic function (p(y,y ,P) on r\A _1 (oo), called the Bloch- 
Floquet function, that has g simple poles Pi(y ), . . . , P g (yo) on T sat- 
isfying X(P k (y )) = >y k (y ). 

This function is a joint eigenfunction of the Schrodinger operator and the 
monodromy operator: 

L s <p(y,y ,P)=\(P)<p(y,y ,P), (2.2) 

<p(y + II, y , P) = fi(PMy, y , P), (2.3) 
and has the following high-energy expansion: 

<p(y, y , P) = exp(t(y - y )VX)(l + o(l)). (2.4) 

The logarithmic derivative of the Bloch-Floquet function is equal to 

, p] = ; <P v (v,Vo,P) = _ i S y (y,X(P)) 

w } ~ <p(y,y ,P) s( y ,\(P)) 2%A(P))' [ - ) 

where the functions R(X) and S(y, A) are defined as 

29 9 
i2(A)=n(A-A,-), S(y,X) = l[(X- lk (y))- (2-6) 

3=0 k=l 

The multi-valued function p(P) = — ^ ln/x(P) is called the quasi-momentum. 
Its differential 

, i da 

dp= -- — 

n \i 

is the unique meromorphic 1-form on T satisfying the following properties: 
1. dp has a single pole of second order at infinity with the principal part 

dp= (-± + 0(l)) dk, 

where k = A^ 1 / 2 is the local parameter, and 
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2. The periods of dp over the a-cycles are equal to zero: 

<f dp = 0, k = l,...,g (2.7) 

Since the function /i(A) = e %Up<yX) is single-valued, the 6-periods of Vt p are 
integral multiples of 27r/IT: 

£dp=^^, n k EZ, k=l,...,g. (2.8) 

Conversely, given a hyperelliptic Riemann surface T together with a two- 
sheeted covering A : T — > CP 1 with real- valued ramification points A , . . . , X 2g 
and oo, and a nonspecial divisor D = P% + • ■ • + P g satisfying A(Pj) G 
[A2J-1, A2j], there exists a smooth real-valued potential u(y) of the Schrodinger 
operator with spectral data {Aj, Pj(0)}. This potential is given in terms of 
the theta-function of T by the Matveev-Its formula: 

u(y) = -2d 2 y \n9iyU! - A(P-y) A(P g ) - K\B i3 ) + C(T), (2.9) 

where A is the Abel map, K is the vector of Riemann constants, C(T) is a 
constant, and the vector U\ is the vector of 6-periods of the unique meromor- 
phic differential dp on Y satisfying properties (1) and (2) above: 

0i)k = ^(f dp k = l,...,g. (2.10) 
2rr Jb k 

This potential is periodic with period II only if the components of the vector 
U\ are integral multiples of 1/11. However, for generic spectral data the 
components of U\ are arbitrary real numbers, so the potential u(y) is in 
general quasi-periodic. 

Therefore, the problem of describing all finite-gap potentials of period 
II is reduced to the following: describe all hyperelliptic Riemann surfaces 
such that the 1-form dp, uniquely determined by conditions (1) and (2), has 
6-periods that are integral multiples of 2tc/H. An effective solution of this 
problem was given in [B], which we now recall. 
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3 Isoperiodic deformations for the Schrodinger 
operator 

Let T be a hyperelliptic Riemann surface corresponding to a real-valued 
potential of the Schrodinger operator of period one. Consider a deformation 
r(t) of T , i.e. let a continuously varying one-periodic family of hyperelliptic 
Riemann surfaces of genus g, equipped with two-sheeted covering maps \(t) : 
T(i) — > CP 1 , branched at points \ (t), . . . , \2 g (t) on the real axis and oo, such 
that T(O) = r . If each of the curves T(t) also corresponds to a potential of 
period one, then this deformation is called isoperiodic. The principal result 
of [B] consists of an effective description of all such deformations. 

Suppose we have an isoperiodic deformation T(t), then each of the curves 
has a unique meromorphic 1-form dp(t) satisfying properties (1) and (2) 
above, such that its 6-periods are integers. Consider the meromorphic 1- 
form 

on r . This form has a double pole at infinity and no other singularities, 
such forms are called weakly meromorphic. If we choose the connection in 
such a way that |j = 0, then the deformation is explicitly given in terms of 
the ramification points 

n\ .,<(\.) 

(3.2) 



dt dp(X k ) 

Conversely, given a weakly meromorphic 1-form ui, we can define a deforma- 
tion of Tq using the above formula. If we now choose the connection in such 
a way that 5| = 0, then we see that 

^ - — (3 3) 

dt ~ d\ [6 ' 6) 

is a single- valued function on r , so therefore the periods of dp(t) are con- 
stant. Therefore, if T corresponds to a potential of period one, then so do 
the surfaces T(t). Thus, isoperiodic deformations of the Schrodinger operator 
are described by meromorphic 1-forms with prescribed singularities, namely 
with a double pole at infinity. We now try to adapt this approach to obtain 
isoperiodic deformations for the acoustic operator. 
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4 Spectral theory of the acoustic operator 

We now consider the acoustic operator 

L A = -r 2 (x)f- 2 , (4.1) 
dx A 

where the function r(x), called the density, is smooth and positive. The 
spectral theory of the acoustic operator with a smooth periodic density has 
been widely studied. The L 2 (R) spectrum of the problem 

L A ?P(x,\) = XiP{x,X), (4.2) 

like that of the periodic Schrodinger operator, has a zone structure and con- 
sists of an infinite number of bands [Ao, Ai], [A2, A3], . . ., with A2j < A^j+i < 
A2J+2. If the density r(x) is smooth, then there is a constraint Ao = 0. The 
finite-gap densities of La can be constructed using the following well-known 
relation [12J: 

Proposition 1 Suppose r(x), ip(x, A) satisfy the acoustic equation 

d 2 

-r 2 (x)—^(x,\) = \iP(x,\). (4.3) 
ax 

Perform a change of variables 

f x dr' 

v(*) = / ^T V (4-4) 







Then the functions 



1 1 

u (v) = j r l( x ) ~ ^r(x)r xx (x), (4.5) 

<p(y,\)=^(x J \)r- 1 / 2 (x), (4.6) 
satisfy the Schrodinger equation 

d 2 

- ^aV(j/) A ) + u (vMy, A) = My, A). (4.7) 

// r(x) is a periodic density of La with period T , then u(y) is a periodic 
potential of L$ with period 

f X0+T dx 

11= / (4.8) 



and the operators La and L$ have the same spectrum in L 
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This proposition allows us to construct finite-gap densities of the acoustic 
operator by taking a finite-gap potential u(y) of the Schrodinger operator, 
given by the Matveev-Its formula for some spectral data {Aj,7fc} with Ao = 
0, and solving the equations (I4.4p - fl4.5l) for the density r(x). This inverse 
transformation is only defined up to a parameter, as the acoustic equation 
has the gauge transformation 

x — > ax, r(x) — > a~ r(ax), (4.9) 

so that a periodic potential corresponds to a one-parameter family of periodic 
densities. The choice of this a should be considered as an additional spectral 
parameter of the problem. 

Explicit formulas for r(x) in terms of theta-functions were obtained by 
Dmitrieva in [2]- [I] by extending the Hopf transformation (14. 41) - 04. 6 p to a 
transformation between the Harry Dym hierarchy and the Korteweg-de Vries 
hierarchy. Finite-gap periodic densities of the acoustic operator are then 
equivalent to x-periodic solutions of the Harry Dym equation 

r t = r 3 r xxx . (4.10) 

Since periodic densities correspond to periodic potentials, and periodic 
potentials of the Schrodinger operator are described by the isoperiodic de- 
formations, the problem of describing all periodic densities of the acoustic 
operator is in some sense solved. However, we would like to find a natural 
choice of the constant for the gauge transformation, depending explicitly on 
the spectral data, and a set of differential equations on the spectral data of 
the acoustic operator, such that the period of the unique density correspond- 
ing to the spectral data is preserved under these flows. We now turn to this 
problem. 

Let r(x) be a finite-gap density of period T, and let i/j±(x, x , A) be the 
Bloch-Floquet function of the corresponding acoustic operator La, that is a 
joint eigenfunction of the acoustic and monodromy operators: 



L A ip±(x, x Q , A) = \ip±{x, x , P), (4.11) 

V±(z + T, x , A) = e ±lTq{x) ^ ± {x, x Q , A), (4.12) 
normalized by the relation 

^ ± (x,x ,T)\ x=xo = l. (4.13) 
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Let L s be the associated Schrodinger operator given by (14. 41) - (14. 6p . Then the 
function ip(x,x , A) can be expressed in terms of the Bloch-Floquet function 
(Q-Q of L s as follows: 

V±(x,x ,A) =r- 1 / 2 (x y/ 2 (x)v(y,y ,HP)). (4.14) 

Hence, we can consider the Bloch-Floquet as a meromorphic function on the 
spectral curve T of Ls, i.e. we can set 

ip±(x, x , A) = ip(x, x , A(P)) (4.15) 

for some meromorphic function ip(x, Xq, P) on T satisfying 

L A iP{x, x , P) = \{P)ip{x, x , P), (4.16) 

ip(x + T, x , P) = e iT(?(P V(z, x , P), (4.17) 

where the multivalued function q(P) is called the quasi-momentum of the 
acoustic operator La- It can be expressed in terms of the logarithmic deriva- 
tive of the Bloch-Floquet function 



as follows: 



1 



xo+T 



g(P) = -/ £(.x,P)dx. (4.19) 

To study the function £(x, P), we first see that it can be expressed in terms 
of the logarithmic derivative x(j/, P) of (p(y,yo,P) as follows: 

ax , P) = (4.20) 
Substituting this in (I2.5p . we get 

?l ' ' r{x)S{x,\) 2r{x) 2 S(x, A) ' 1 ' } 

29 g 

where P(A) = A J^[(A — \j) and S(x, A) = J^[(A — 7fc(y(x))). Therefore, 

j=i fc=i 
the quasi-momentum q(P) and its differential have the following high-energy 
expressions: 

^) = ^f +T |) +0(1) - (422) 

9 



/ i r o+T dx \ 
d " {P)= {-¥fl W) + ° m ) Ak > (423) 

where k = A -1 / 2 . 

A simple calculation also shows that £ (x, P) satisfies the Ricatti equation: 

- l ^'( x ,P)+e(x,P)-^p- = 0. (4.24) 

The spectrum of the acoustic operator always starts at A = 0, i.e. the left- 
most branch point of V is A = 0. We consider the Ricatti equation in the 
neighborhood of this point. Let k = y/\ be the local parameter and consider 
the Taylor series for £(x, P) at k = 0: 

C(x, P) = £ (x) + &{x)k + U*)k 2 + 0(k 3 ). (4.25) 

Substituting this into the Ricatti equation, we get that £,o(x) = 0, £i(x) = C 
for some constant C, and that £2^) is a full derivative. Therefore, the quasi- 
momentum q(P) and its differential near k — are equal to 

q(P) = Ck + 0{k 3 ). (4.26) 

dq(P) = (C + 0(k 2 ))dk. (4.27) 

On the other hand, comparing the obtained expression for P) with (I4.2ip . 
we see that 

C = 7 r, ^ (4-28) 

and hence r(x) is expressed in terms of the spectral data and the additional 
constant C as follows: 



r(x) = , ^ Xl '/' X29 n . (4.29) 

The constant C in this formula corresponds to the gauge transformation (14.91) 
and should be considered in addition to the Xj and 7 ^ as part of the spectral 
data of the problem. 

To construct isoperiodic deformations of the acoustic operator, we need a 
rule of choosing the constant C. The natural choice seems to be C = (— l) 9 . 
With this choice, the value of dq at A = does not depend on the spectral 
curve T, which will be used to construct the isoperiodic deformations. Also, 
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with this choice of C, a vanishingly small potential u(x) — > 0, or equivalently 
vanishingly small energy gaps (A 2j -i — A 2j - — > 0), correspond to a density 
r(x) — > 1 normalized at unity. 

We summarize the results of this section in the following 

Theorem 1 Let {Aj, 7fc(y)} fre £/ie spectral data of a finite-gap periodic Schrddinger 
operator Ls with potential u(y) such that Ao = 0. Then the density r(x) of the 
associated acoustic operator La related to u(y) by ( |^.^[ )-( f^T^ is given by the 
equation h4-2ty , where C is an arbitrary constant. If we choose C = (— I) 9 , 
then the differential of the quasi-momentum of La has a pole of second order 
at A = oo with principal part ffj^fy , and has fixed first and second order 



terms at A = 0: 

dq(P) = ((-iy + 0(k 2 ))dk. (4.30) 

Using this theorem, we construct isoperiodic deformations for the acoustic 
operator. 



5 Isoperiodic deformations of the acoustic op- 
erator 

We recall that a deformation of a Riemann surface is given by a meromorphic 
1-form 

LU = ^d\. (5.1) 

dt v ; 

To give an isoperiodic deformation for the acoustic operator, this 1-form must 
satisfy the following conditions: 

1. The quasimomentum differential dq has a second order pole at infinity, 
so || has a first order pole at infinity. Since d\ has a third order pole 
at infinity, u has a pole of fourth order at infinity. 

2. At zero dq has fixed first and second order terms, therefore || has a 
third order zero (with our choice of C). Since dX has a first order zero, 
uj has a zero of fourth order at A = 0. 



Therefore, a 1-form u defines an isoperiodic deformation if and only if it 
has divisor 4 • — 4 • oo. The space of such forms is (/-dimensional. Since the 
total space of spectral curves has dimension 2g, and there are g conditions 
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for the period to be equal to one, these forms are a basis of isoperiodic 
deformations. 

As before, the deformation is given explicitly in terms of the branch points 

as 

j = l,...,2g. (5.2) 



dt dq(Xj) 

On a hyperelliptic surface, the 1-form dq can be written explicitly as 

d ^T7m dx ' (5 - 3) 

where 

Q(X) = qg X3 + ... + q , (5.4) 

and q = (— 1) 9 WX[ ■ ■ ■ X 2g - The coefficients of the polynomial Q(X) are 
determined by setting the a-periods of dq to zero and are expressed in terms 
of certain hyperelliptic integrals. 

An arbitrary 1-form uj with divisor 4 ■ — 4 • oo can be written as 

u = f i-)— dX, (5.5) 

where 

/(A) = f g+1 X° +1 + ... + f 2 X 2 = X 2 (f g+1 X^ + ■ ■ ■ + f 2 ). (5.6) 
The deformation given by such an u has the form 
dX l __f{X^ 

dt ~ Q(X 3 y 3-^---^9- (5.7) 

It is possible to choose a basis in the space of meromorphic 1-forms with 
divisor 4 • — 4 • oo and construct a basis of deformations by this formula. 
However, the right-hand side of the deformation equations contains the coeffi- 
cients of Q(X) that are in turn expressed as hyperelliptic integrals containing 
the ramification points. To avoid this difficulty, we use the approach of [6]. 
Factoring the polynomial Q(X) 

Q(X) = {faX - y^X 2 ) ■ ■ ■ {f3 g X - v^-^), (5.8) 

we consider the parameters fa as independent variables that are deformed 
along with the ramification points. This extension of the phase space greatly 
simplifies the deformation equations. 

The deformations of fa are easily shown to be given by the following 
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Lemma 1 Let be the flow given by the form u with polynomial /(A) . Then 
the deformations of Aj and fl k have the form 



dt 



= I3 k 



dt Q(\j) 
1 /(A 2fc _i) 1 



(5.9) 



/(A 



2k j 



2A 2 fc_i <5(A 2 jk_i) 2A 2 fe Q(A 2 fc) 
/(A)#(A)' 



Q'(A) 



(/'(A) - 



2i2(A) 



■ l _V /A 2fc-l A 2fc 



(5.10) 



As we see, the equations no longer contain hyperelliptic integrals. This 
form makes the equations of isoperiodic deformation useful for numerical 
simulations. 

To write the deformations explicitly, we choose a basis in the space of 
differentials with divisor 4 • — 4 • oo as follows: 

CfcA 2 



AcA — A/A 2 fe_iA 2/ t 
where Ck are arbitrary constants. Then 

/fc(A) c fc A 2 



dq, 



and the deformation equations are given by our final theorem. 



(5.11) 



(5.12) 



Theorem 2 Let T be a hyperelliptic Riemann surface with ramification points 
A = 0, Ai, . . . , A 2g; corresponding to a smooth density of period one. Let the 

zeroes of the quasimomentum be — ' ^p^* - , k = 1, . . . , g. Consider the flows 

(5.13) 



dt k 



I3 k \j — a/ A 2 fc-iA 2 fc 



w 

dt k 



77- = C k f3 k 



A 2 



21-1 



\2 
A 2l 



2A a _i (3 k \ 2 i-i - A/A 2fe _iA 2 fc 2A 2 / f3 k \ 2l - a/A^-iA^ 



A 2 ;-iA 



21 



01 



a/ A 2 z-iA 

A 



21 



— a/ A 2 fc-iA 



2fc 



I ± k 



(5.14) 



13 



dt k 



\ 2 

A 2k-1 



\2 

A 2k 



2A2fc-l (3k\2k-l — a/ ^2k-1^2k ^21 f3 k X 2 k — a/ A2fc-lA2fc 



+ 



1 I A2/t-iA 



2k 



l^k 
A2/C-1A2A; 



$> [Pi 

l+k \ 



2k-l^2k 



0k 



\/ A2«-lA 



21 



a/ A2A:-lA2fc 



29 

E 



2fc-lA2fc , =1 a/ A2fc-lA2fe — Xjftk 



(5.15) 



on £/ie spectral data, where ^- is the flow given by Uk defined by A5.11\) . 
Then any Riemann surface obtained by moving along these flows in the mod- 
uli space of hyperelliptic Riemann surfaces corresponds to a density r(x) of 
period one of the acoustic operator. These flows form a basis in the space of 
isoperiodic deformations of the acoustic operator. 

We see that these extended equations do not involve hyperelliptic inte- 
grals. 
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